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'^ . Abstract 

The average charged hadron multiplicity in the e~^e~ events with 
C^ ■ the primary it-pair at the coUision energy 500 GeV, as well as the aver- 

J^ . age multiplicity of charged hadrons from the top quark are calculated 

^ ; in QCD to be 86.7 ± 1.11 and 41.0 ± 0.54, respectively. 

cn _ 

1 Introduction 

!>■ . 

O . Experiments at LEP and SLAG revealed, besides other important results, 



quite interesting feature of the hadron multiple production dependent on 
the mass of the "primary" (anti) quarks which launch the process of the 



^ ■ QCD evolution. It appeared that differences between the light and heavy 



quark- induced multiplicities become energy- independent. QCD calculations 
describe the phenomenon quite well. 

Certainly, LEP could not give the information on the events induced by 
the top quarks. Recent discussions of the ILC project give us occasion to 
provide QCD predictions concerning the hadron multiple production in the 
events with primary t-quarks. 
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We manage to calculate the average hadron multiplicity Nti in the e^e~ 
events with ti pair at the collision energy of the ILC with the following 
prediction: 

Nu{W = 500 GeV) = 86.67 ± 0.55 . (1) 

We also theoretically calculated the average hadronic multiplicity from the 
top quark: 

nt = 41.03 ±0.27. (2) 

Both values correspond to the average value of the top mass rrit = 170.9. 
Everywhere below, it is assumed that we deal with average multiplicities of 
charged hadrons. 

The paper is organized as follows. In order to make our calculations of 
the hadron multiplicity in top quark events more easy for understanding, we 
consider first the multiple hadron production in cc (bb) events. The hadron 
multiplicities in e~^e~ events associated with the tt-pair production are cal- 
culated in Section [3] in the framework of perturbative QCD. In Section H] the 
numerical estimations and our main results are presented. 

2 Hadron multiplicity in e+e~ annihilation 
associated with cc or bb-pair production 

Hadron multiplicity in qq event, Ngg{W), can be represented in the following 
general form [Ij: 

Ng,{W) = 2n, + C,l^^n,{e)Eg[^) , (3) 



"o 



where q means a type of quarks produced in the process of e~^e~ annihilation 
into hadrons at the collision energy W. In what follows, the notation q = Q 
{heavy quark) will mean charm or beauty quark, while the notation q = I 
[light quark) will correspond to a massless case (when a pair of u, d or s- 
quarks is produced, whose masses are assumed to be equal to zero). The top 
quark production (g = t) will be studied in Sections |3] and HI 

The first term in the r.h.s. of Eq. ([3]), 2rag, is the multiplicity of primary 
(anti)quark of the type q (i.e. the multiplicity from the leading hadron which 



contains this (anti)quark). It is taken from an analysis of the data (2nc = 5.2, 
2nfe = 11.1 [3J, and 2n, = 2.4 [3]). 

The quantity ng{k'^) in ([3]) is the mean multiphcity of the gluon jet with 
a virtuahty /c^, for which we will take a QCD-based parametric form, with 
parameters fit to data, while Eq{k'^/W'^) is the inclusive spectrum of the gluon 
jet emitted by primary quarks. It was explained in detail in Ref. [1] that one 
should not consider this mechanism of hadron production via gluon jets as 
due to "a single cascading gluon". That quantity E{k'^/W'^) is an inclusive 
spectrum of the gluon jets is seen, e.g., from the fact that the average number 
oiiets J dk^/k^Eq{e/W^) ^ 1. 

Qo is a phenomenological parameter denoting the scale at which "precon- 
finement" of the off-shell partons occurs (as explained in Ref. [S]). 

Let us introduce variable^!! 

7? = In— (4) 



and 



as well as notation 



Ql 



y = ^^T^^ (5) 



n 



Cpttsik^ 



'~ TT 



■n,(fc2), (6) 



where Cp = {N^ — l)/2Nc, and A'"c = 3 is a number of colors. Then Eq. 
can be represented as 



Y 



Nqq{Y) = 2nq+ fdr] hg{Y - r/) Eq{r]) =2nq + iV,(F) . (7) 



In particular, Nii{Y) means the multiplicity of hadrons in light quark events, 
while Nqq{Y) denotes the multiplicity of hadrons in a process when a pair 
of the heavy quarks is produced. 

The physical meaning of the function 

Y Y 

Nq{Y) = fdr]ng{Y - r])Eq{7]) = /"rfr/'n,(r/') Eq{Y - r]') (8) 





^We will often use "rapidity-like" variables (analogous to rj and Y) instead of the energy 
variable W throughout the paper. 



is the following. It describes the average number of hadrons produced in 
virtual gluon jets emitted by the primary quark and antiquark of the type 
q. In other words, it is the multiplicity in qq event except for multiplicity of 
the decay products of the primary quarks at the final stage of hadronization 
(the terms 2nq in ([7])). 

For the massless case, the function E = Ei was calculated in our paper [T] . 
In terms of variable 

a = exp{-r]) , (9) 

it looks as 



E[r^ia)] 



(1 + 2cr + 2a^) In - 



1 3 + 7a 



cr 



ail + a] 



a 



a 



TT 



+ 4(t(1 + a) — + In a ln(l + a) + Lial-a) 



(10) 



where Li2(-2) is the Euler dilogarithm. The function E{r]) is presented in 
Fig. [H It has the asymptotics 
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Figure 1: The function E[y). 



The derivative of E{ri) is positive, and dE{r])/dr] = at ?7 = 0. As 
a result, the associated multiplicity Nq(W) ([H]) is a monotonic increasing 



function of the energy W for any positive function ng{k'^) since 

dNAY) /"^ , , , dE(Y -7]) , , 

Now consider the multiphcity difference in events with the hght and heavy 
flavors (Q = c or 6): 

^Qi = NQQ-Nn. (13) 

At ly ^ m-Q, one can neglect small power-like corrections Qirn? jW'^\ In 
such a case, the quantity bqi is defined by [1] 

5Q, = 2(nQ-n,)-AiVQ(FQ), (14) 

where the notation 

^Nq{Yq) =Ni-Nq= jdyh.iYQ - y) AEgiy) , (15) 



as well as variables 



and 



2 



?/ = ln^ (16) 



2 

Yq=^^J^ (17) 



are introduced. The lower limit of integration in Eq. flTSl) . — ln(iy^/mQ), is 
taken — cxd because of the fast convergence of the integral at negative y. 
Let us use another dimensionless variable 

p = exp(-y) . (18) 

The explicit form of AEq was derived in our paper [T] (see Fig. [2]): 

AEQ[y{p)] = [l + p (Ip - 3) In 1 + (^ + 7p) 

+ p (7p - 20) J(p) + 20 i^^] . (19) 



where 

J{p) = < 1 , P = 4 , 

The function AEq^y) decreases at y ^ — cxo (p -^ oo) as 

y—>—oo O 

and has the following asymptotics at y ^ oo (p ^ 0): 



Ai?Q(y) 



1 

y- o • 



(20) 
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(22) 



y^co 



AEQ(y) 




Figure 2: The function AEq^y). 

Thus, we get the relation between average multiplicities of hadrons in QQ 
and // events: 



NggiW) = NaiW) - ANQimg) + 2{nQ - m] 



(23) 



with the multiplicity difference ANq defined by Eq. flTSj) . 

Our calculations [Ij of the multiplicity differences 6qi = Nif — Nqq {Q = 
b, c) with the use of formula ( 123|) appeared to be in a good agreement with 



the data. Recently we have reconsidered the QCD upper hmit on quantity 
5uf^ which appeared to be very close to all present experimental data on 

3 Hadron multiplicity in e+e~ annihilation 
associated with tt-pair production 

The goal of this paper is to calculate A^^j, the average multiplicity of hadrons 
produced in e^e~ events with the primary ti pair. We consider the case 
when the top (antitop) decay mode is pure hadronic. As a byproduct, we 
will calculate nj, the hadron multiplicity of the on-shell top decay products. 
We will assume that the square of the matrix element of the process 
e^e~ — >■ t*i* — > X is factorized as follows: 

|M(e+e" -^ ft* -^ hadrons) p = |M(e+e" -^ ft* ^ ti + hadrons) 1^ 

X \M{t^ hadrons) p 
X \M{i^ hadrons) p , (24) 

where t*{t*) denotes the virtual top quark (antiquark). 

The factorization of the matrix element fl24p means that there is no signif- 
icant space-time overlap in the decay products of the on-shell t and t-quarkso 
The QCD non-singlet evolution of the primary virtual t-quark is very slow 
because the difference of virtualities in logarithmic scale is very small down 
to the top quark mass. In other words, the virtual t-quark becomes "real" 
after just a few gluon radiation. 

The effect of possible color reconnection was investigated by comparing 
hadronic multiplicities in e^e~ — >■ W^W~ — >■ qq'qq' and e^e~ -^ W^W~ — >■ 
qq'lvi events. No evidence for final state interactions was found by measuring 
the difference (^4^) — 2(n2g;p) [H [8]. From the space-time point of view W 
bosons and t-quarks behave in a similar way, i.e. the latter manage to cover 
the distance M ~ l/Ff, where Vt is the full width of the top. Since Vt ~ Vy/i 
we expect no interference effects in the decays of the on-shell t and t-quarks. 



^The data on Na as well as on Su at different energies corrected for detector effects as 
well as for initial state radiation were recently cited in [5]. 

•^Note that the off-shell t and i-quarks fragment into hadrons through the emission of 
the gluon jets in a coherent way (the first term in the r.h.s of Eq. 



According to Eq. 
the formula: 



the associative multiphcity in tt-event is given by 



where 



N^t{W,mt)=Nt{W,mt) + 2nt, 



{W-2m.tf 



N,{W.^t) = CF I ^^^n,{e)E,{W',k 



2 2\ 



(25) 



(26) 



Qg 



Here and in what follows we will assume that the collision energy is a 
typical ILC energy, W = 500 GeV, for definiteness. In such a case, contrary 
to Eq. (IH]), power corrections Olrrit/W) should be taken into account. The 
explicit form of the inclusive distribution of the gluon jets with the invariant 
mass V^ looks like 



Et{q',e,m^ 



"dk 



+ 2ri- 



In 
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'l(g2 + p)2_4^4 


o^ 


J] g4 
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v + V 


^2-l^(A-r^)/(Ao 


-v) 



k q^ + k'^ + 2m2 



r^- v^r72-l^(A-r/)/(Ao-r7) 



X 



(l + 2mVg2)(l + 2mVA;2) 



q- 

1 + 



where k 



k"^, q 



A 



/q"^, and the following notations are introduced: 
q^ + k"^ - Am^ q^ + k'^ 



(27) 



2qk 



An 



2qk 



(28) 



This formula has been derived by calculating QCD diagrams in the first 
order in the strong coupling constant jj In the massless case {m = 0), 
we immediately come to the function E{k'^q^) flTO|) . while by neglecting 
small corrections 0(m^/g^), one can derivqj the explicit form of the func- 
tion ^EQ{e/m^) = EQ{q^,k^,m^) - E{q^,k'^) 1^. In our case (g^ = VV^, 



^See our comments after Eq. ([3|). 
""After variables are properly changed. 



m = nit) "we will estimate the integrals in Eqs. (1261) . (1271) numerically (for 
details, see Section Hj). 

Now let us calculate another quantity in Eq. fl25|) . rit, which describes the 
hadronic multiplicity of the t-quark decay products. The top weakly decays 
into W^ boson and 6-quark. In its turn, the W^ boson decays into a quark- 
antiquark pairjj The quark-antiquark system results in massive jets which 
fragment into hadrons (see Fig. [3]). 




+ 




Figure 3: The emission of the gluon jets (spiral lines) by the quark pair 
resulting from the decay of the VT^-boson in the first order in the strong 
coupling constant. The W~^ boson is produced in the weak decay of the top. 



The gluon jets can be also emitted by the on-shell t-quark before its 
weak decay (the first diagram in Fig. Hj) or by off-shell bottom quark (the 
second diagram in Fig.H]). At the end of these emissions, the on-shell 6-quark 




-t- 




Figure 4: The emission of the gluon jets by the on-shell top quark and by 
off-shell bottom quark in the first order in the strong coupling constant. The 
emissions take place before and after the weak decay of the top, respectively. 
Off-shell quarks are denoted as t* and b*. 

weakly decays into hadrons whose average multiplicity is equal to nf,. Since 



^Remember that we are interested in hadronic decays of the W boson. 



VT-boson is a colorless particle, the diagrams in Fig. H] do not interfere with 
those presented in Fig. [31 

Thus, the multiplicity nt is a sum of three terms: 

nt = nw + nth + n^ . (29) 

The quantity n^ is experimentally measurable one [3]. The first term in 
Eq. (!29l) . niy, is the hadron multiplicity of the W boson decay products. 
The second term, ritb-, is the hadron multiplicity in the gluon jets emitted by 
the on-shell top quark before its weak decay as well as by the bottom quark 
after the top decay. 

3.1 Multiplicity of W boson decay products 

The W^ boson can decay either into two light quarks {ud and us pairs) or 
into cd (cs) pair. The former case is treated analogously to the light quark 
event in e+e~ annihilation taken at the collision energy W = mw- Here we 
will study the latter case. 

Let Nqi{W) be hadronic multiplicity associated with the production of 
one heavy quark (antiquark) of the type Q and one light antiquark (quark) 
of the type /: 

NQiiW) = (UQ + m) + NgiiW) . (30) 

Now let us introduce the notation (not to confuse with ANq from above): 

ANqi = Ni- Nqi . (31) 

Then the first term in the r.h.s. of Eq. fl29l) is given in terms of the function 
ANci by the formula: 

nw = NiiiYw) + - [-ANa{Y,) + n, - m] , (32) 



where Y^ = ln{ml/Ql) and 



Yw = ln^. (33) 

Qo 



The function Nii{Y) in fl32l) is the hadronic multiplicity in light quark events. 

Thus, we need to find an expression for ANd- Note that the formulae 

( TT5|) . ( TT9l) from Section [2] correspond to the case when a pair of heavy or pair 

10 



of light quarks is produced. Now we have to study the case when hadrons 
are produced in association with a single heavy quark (namely, c-quark) and 
one light quark. 

Our QCD calculations result in the following representation for the mul- 
tiplicity difference (see Appendix for details): 



ANqi{Yq) = / dyngiYg - y) AEgiiy) 



(34) 



with the dimensionless function AEqi: 



AEQi[y{p)] = i [2 + p (3p - 2)] In - + i (5 + 6p) 
+ ip(3p-8)J(p) + 6i-^. 



(35) 



Here p = exp(— y). The quantity J{p) was defined above fl20l) . The function 
AEqi^y) is shown in Fig. O 

AEQi(y) 




Figure 5: The function AEQi{y). 



Since 



AEQiiy) 



y—*—oo 2 



(36) 
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the integral fl34p converges rapidly at the lower limit. Asymptotics of AEQi{y) 
at large y is the following: 



^Egiiy) 



y-*oo 



2 1^-2 



(37) 



AEq(y), 2AEqi{y) 




8 10 



Figure 6: The function AEq^y) (solid line) vs. function 2AEQi{y) (dashed 
line) . 

We derive from Eqs. ([22]), ([37]) that AEgiiy) = 0.5AEq(|/) at large y. 
Numerical calculations show that 2AEqi is very close to AEq at all y (see 
Fig. [6]). Thus, we can put for our further numerical estimates: 



AN,i = IaN,. 



(38) 



This relation means that 



NqI = I [Nu + Nqq] = Nij +15qi. 



(39) 



Correspondingly, we obtain: 



nw = Nfiimw) + t 4; 



(40) 
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3.2 Multiplicity of top and bottom decay products 

As was already said above, the on-shell top quark can emit jets before it 
weakly decays into W'^b. After the weak decay of the top, the off-shell 
6-quark "throws off" its virtuality by emitting massive gluon jets. The frag- 
mentation of these massive gluon jets into hadrons results in the average 
hadron multiplicity rifb- 

To calculate the multiplicity ntb, one has to derive the inclusive spectrum 
of the gluon jets, emitted by the top and bottom quarks. Let us denote it as 
Etb- Then the multiplicity ntb will be given by the formula: 

Ytb 

ntb = / dyfigiYtb - y) Etb{y) , (41) 



where 

y ^ In ''"'-';»'-'"''' ■ (12) 

and 

{mt -mw- rubf , . 

Yth = In -r^ , (43) 

with k"^ being the gluon jet invariant mass, {mt — myy — mb)'^ its upper bound. 
In the lowest order in the strong coupling constant, the quantity Etb{y) 
is given by two diagrams in Fig. HI It is presented by an integral which 
depends on the ratio /c^/m^, as well as on mass ratios m^/ml and ml/m^. 
This integral cannot be calculated analytically, but can be estimated numer- 
ically. The function Etb{y) is presented in Fig. [71 It is worth to note that in 
the Feynman gauge the dominating contribution to Etb{y) comes from the 
interference of two diagrams shown in Fig. |H 

3.3 Associated multiplicity of hadrons in tt events 

The formulae of the previous subsections enable us to derive the average mul- 
tiplicity of the charged hadrons in e~^e~ annihilation at the collision energy 
W associated with the production of the tt-pair. It is of the form: 

N,){W, mt) = Nt{W, mt) + 2[Na{mw) + ntb + n^] 

+ [-/\Nd{mc)+nc-ni]. (44) 
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Figure 7: The function Etbiy). 

Let us remind to the reader the meaning of all quantities in Eq. (144 
The function Nt{W,mt) describes the average number of hadrons produced 
in association with the tt-primary pair, except for the decay products of the 
top(antitop) fl26|) . The quantity Nijijnw) is the mean hadron multiplicity in 
the light quark event taken at the energy E = mw- The hadron multiplicity 
nth comes from the emission by t and b quarks fj4Tl) . The quantity ni is the 
mean multiplicity of hadrons produced in the decay of the on-shell primary 
quark q [q = I, c, b). Finally, the combination [ANci{mc) + rii — ric] is the 
difference of multiplicities in the processes with the primary //- and c/-pairs. 
As for the hadron multiplicity resulting from the decay of the on-shell top 
(anti)quark, it is given by 

nt = nw + ritb + n^ 

= Nfi{mw) +nth + nh + - [-AA^c/("^c) + r^c - rii] . (45) 

The expressions for iV;;-, AiVc^ are given by Eqs. ([7]), (IMI) . respectively. 

Let us stress that Nij and Uq (g = 6, c, /) are extracted from the data, and 
AA^cZ can be related with the measurable quantities (see formulae (l38l) - (H0ll 
in the end of subsection 13.11) : 

-AN, = n,-ni-- 6ci . (46) 
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Then we obtain: 



and 



nt = Niiimw) + T 4/ + ntb + n^ , 
<-(iy, mt) = Nt{W, mt) + 2 \Nfi{mw) + ^^01 + ntb + rn, 



(47) 



(48) 



where ntb is defined above fj4Tl) . In what follows, we will use the value 

5d = 1.03 ± 0.34 (49) 

from Ref. [6]. 



4 Numerical estimates of hadron multiplici- 
ties. 

In order to estimate the multiplicity of the decay products of the top (for- 
mula (H71) ). one has to know the energy dependence of the hadron multiplicity 
in the light quark event. The latter is defined by Eq. (JTj), where the function 
fbg is related with gluon jet multiplicity Ugiji^) ([6]). We have fitted the data 
on NiiiyV) by using the following QCD-motivated expression for ng{k'^): 



Ugik ) = a + 6exp cJ\\i{k'^/Ql 



(50) 



where c = 1.63, and k"^ is the invariant mass of the jet. We have got the 
following values of the parameters: 



3.89, 



h = 0.01, Qo = 0.87 GeV . 



(51) 



The result of our fit is presented in Fig. [H] in comparison with the data. 
Note that y^ /d.o.f. becomes twice smaller if one eliminates the experimental 
point at FT = 58 GeV (open circle in Fig. [8]), which lies much lower than 
neighboring points, is crossed out from the fit. In such a case, the values of 
the parameters are practically the same as in (IST!) with x^/c/.o./. = 0.90. 

For our numerical estimates we shall use my/ = 80.40 ± 0.03 GeV [lU] 
and the recent value of the top mass [H] : 



TUt 



170.9 ±1.8 GeV. 



(52) 
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Figure 8: The fit of tlie data on tfie liadron multiplicity in liglit quark events 
(solid line, see formula in the text). The data are taken from Ref. [6]. 

As for the bottom quark, its pole mass is quoted in [TU] to be m^ = 4.7 — 5.0 
GeV. We will take the average value 



mb = 4.85 ±0.15 GeV. 



(53) 



By using our fit, we obtain Nij^mw) = 19.09. Then we get (see Eqs. (j^Oj) . 
dlHD, gl])): 

nvK = 19.34 ±0.10, (54) 

nt6 = 16.14 ±0.24. (55) 

The error in Eq. (!54|) is defined by that of the multiplicity difference 6ci (149|) , 
while that in Eq. (I55p comes from uncertainties of the quark masses rrit (I52j) 
and rrih ( 155]) . 

Our result flMl) is in a nice agreement with the experimental values from 
Ref. |7|, 

(56) 



and Ref. [8], 



nw = 19.3 ±0.3 ±0.3 



nw = 19.44 ±0.13 ±0.12 



(57) 
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Now let us calculate the associated liadron multiplicity in ti event (148 p at 
fixed energy W = 500 GeV. To do this, we need to estimate the multiplicity 
Nt{W,mt) by using formulae fl26|) and ([2] 



Nt{W = 500 GeV) = 4.61 ± 0.11 . (58) 

The errors in fISS]) come from top quark mass errors. It follows From Eqs. flM]) . 
dSHD that Ntiie+e" ^ ti ^ hadrons) = 86.67 ± 0.55. 

In order to estimate possible theoretical uncertainties, we have repeated 
our calculations taking the different form of the average multiplicity of the 
gluon jet (compare with the QCD-based expression (150|) ): 

hg{e)=A + B\n^^. (59) 

It appeared that the data on the average multiplicity in light quark events 
can be fitted well by using this expression (with A = 4.21, B = 0.012 and 
Qo = 0.93 GeV). In particular, we have obtained the following average values 
for the hadronic multiplicities: nw = 19.52, nth = 16.43, Nf = 4.59. Thus, 
theoretical uncertainties can be estimated to be 0.47 and 0.96 for rit and N^t, 
respectively. 

Taking into account the phenomenological value of n^ [6] , 

rib = 5.55 ±0.09, (60) 

we obtain from fj47|) : 

nt{t -^ hadrons) = 41.03 ± 0.54 . (61) 

In the case when the W boson decays into leptons, we predict: 

nt{t -^ Wi + hadrons) = 21.69 ± 0.53 . (62) 

As a result, we obtain the average hadron multiplicity in e~^e~ annihilation 
with the primary tt-pair: 

Nti{e+e- ^ti^ hadrons) = 86.67 ± 1.11 . (63) 

In the case when both W bosons decay into leptons, the mean multiplicity 
of the hadrons is expected to be 

iVtt-(e+e- ^ti^ W+W- + hadrons) = 47.99 ± 0.59 . (64) 
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Finally, we predict that 

Nttie+e- ^ti^bb W^W + hadrons) = 36.89 ± 0.56 . (65) 

All estimations fl63l) - fl65l) correspond to the collision energy W = 500 GeV. 
We can mention the estimation of the hadron multiplicity from Ref. [H], 
Nu[e'^e~ -^ ti ^ bb W~^W~ + hadrons) ~ 29, which was obtained for 
W = 390 GeV and rrit = 175 GeV. For the same values of W and mt, our 
formulae give Nti{e^e' ^ ti —>■ bb W~^W~ + hadrons) = 34.3. 

The formulae flM]) - fp^ is our main result. We hope that the hadron 
multiplicities of the top decay products (Eqs. flMl) and fl^ ) will be measured 
at the LHC. 

Acknowledgements 

We are thankful to the referee for his comments and critical remarks that 
helped us to improve the presentation of some our results. 

Appendix A 

Here we present some formulae for the case when the W boson decays into 
hadrons via production of cs (or cd) pair. Since the total width of the W 
boson, Tw, is much less than its mass, and the hadron multiplicity is a 
smooth function of energy, we will use zero width approximation and take 
the multiplicity a.t W = m\y. It can be shown that the account of the W 
boson width results in corrections which are numerically small (less than 
1.6%, see Appendix B). 

To calculate the inclusive spectrum of the gluon jets emitted by the de- 
cay products of the PF-boson, we need to calculate two sub-diagrams of the 
diagram presented in Fig. O 

The exterior part of the diagram in Fig. M describes the emission of the 
W boson by the top quark, with / being a 4-momentum of the t-quark, q is 
a 4-momentum of the W boson, and (/ — g) is a 4-momentum of the 6-quark. 
The corresponding expression for this part of the diagram (after convolution 
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Figure 9: The generalized diagram describing the inclusive spectrum of the 
gluon jets (curly line) with the virtuality k"^ inside the W boson (wavy line). 
In its turn, the W boson is a product of the weak decay of the top quark 
(solid line with 4- momentum /). 



with the tensor parts of the W boson propagators) looks like 
'^.^^,y{l, q) = 4:1 41^1^ - 2 {lf,q„ + q^Q - g^,u{ml + ml - 



T^Wj 



m 



[{ml - ml)l^ - mlq^] q^ 
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% [(m^ - ml)l^ - mlq^] 
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+ 
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■ qf^q,^ [( 



TTl, 



m 



2\2 



q\ml + ml)]+2ie^,pJPq'' 



w 



(A.l) 



where mt and mf, are masses of the top and beauty quark, respectively. In 
what follows, we will neglect power corrections of the type Olrric/mt) and 
Oimb/rrit). 

The inner part of the diagram in Fig. [9] describes the distribution of the 
massive gluon jet with the invariant mass k"^ produced by the W boson. Let 
D^'^ be the expression corresponding to this diagram. In the first order in 
the strong coupling constant, D'^'^ is represented by the sum of three QCD 
diagrams presented in Figs. llOal llObI and llOcI 

Since q^H^u = 0, one needs to calculate only two tensor structures in 
Dfj^^, namely, g^^ and kfj_k^. The convolution of the tensor 11^,^ with the 
tensor D^" = g^^Di + fc'^fc'^Ds + ■ ■ ■ , 



A = li^^D^"' = g^-" U^,D,{e, qk) + k^k^ ^^./^2(A:^ qk) 



(A.2) 
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Figure 10a: The inclusive distribution of the massive gluon jet with the 
virtuahty fc^. The wavy hne is the W boson, whose 4-niomentuni is q. The 
thick quark hne is a heavy quark, while the thin line is a light quark. The 
cut quark lines mean that these quarks are on-shell quarks. 



depends on Lorentz-invariant variables k'^,qk, Ik. Moreover, it is a polyno- 
mial of the second order in variable Ik. One can use the following useful 
relation: 



d'^k 
where 



A{k'^,qk,lk) 



(27r)2m2(l-r) 




dPd{qk)d{lk)A{e,qk,lk) 



m 



w 



mt 



(A.3) 
(A.4) 



It is naturally to integrate the function v4(fc^, gfc, Ik) first in variable Ik, 
whose lower and upper limits are 



{lk)^ = ^^\{qk){l + r)±{l-T)\l{qkY 



m'^k^ 



(A.5) 
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Figure 10b: The same as in Figs. I10a[ but with the gluon jet emitted by the 
hght quark. 




Figure 10c: The interference diagram which also contributes to the inclusive 
distribution of the gluon jets with the virtuality k"^ inside the W boson. The 
diagram are taken in the sum of the diagrams with the factor 2. 



by using the following formulae: 

I ilk) 



r{^i^)+ 1 I 

J{ik)- I' ^ 

riik)+ 1 / 

/ (Ik) d{lk) = — (1 - r') (qk) J{qky - ml^k 

J(ik). '^f ^ 



m 
m 

(Ik). 



{ikfd{ik) 



(l-r) A{qkY{l + r + r^)-m^k\l-r) 



12r3 
X J {qky — m'^fk'^ . 



(A.6) 



21 



Note that the parts of tensors Il^,y and D^'^ antisymmetric in indiceqil give 
no contribution after integration in Ik. 

Integration hmits in variable qk looks like 



{qk)+ = I (ml, + e) . (A.7) 

As a result, we obtain the hadron multiplicity associated with the charm 
quark in a hadronic decay of the W boson (with the multiplicity of primary 
quark decay products subtracted): 



-^w iQk)^ 



^^' = i2nrmlil -r)H /^^' [^"^^^'^J / ^^^^^ 

Qg (gfc)- 

m+ 

X f dilk)A{e,qk,lk) . (A.8) 

m- 

The dimensionless quantity ng{k'^) in (1A.8I) is the multiplicity of hadrons in 
the gluon jet whose virtuality is k'^, while the normalization H is given by 

The analytic expressions for the functions Dilk"^, qk), D2{k^, qk) in (]A.2p 
are rather complicated to be shown hereo That is why we present only the 
final results of our QCD calculations based on the formulae of this Appendix 
(see Eqs. ( 1M1) . (1351) in the main text). 



Appendix B 

In this Appendix we will demonstrate that the account of the W boson width 
results in only small corrections to the hadronic multiplicities. 



'^As the last term in Eq. (jA.ip . 

®They depend also on the masses rric and mw- 
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The denominator of the W boson propagator (see the diagram in Fig. ^ 

q^ -m^^ + imwTw , (B.l) 



where rriw is the mass of the W boson, T^i its full width. The mean multi- 
plicity is given by the formula 



where Uh {rih = riw or ritb, see the main text) depends on the W boson 
virtuality q^. Here 

F{q^) = {ml-q^)\ml + 2q^), (B.3) 

and the normalization A^ is 

oo 



2 TJlw^W I?,'^2^ 

(g2 - m^)2 + ml^Tl. 





N- d,^ ,, 'T" ... F{f) . (B,4) 



In both (n/i) ( ]B.2p and A^ (IB. 31) the factor niw^w is introduced, while com- 



mon constants are omitted. 
In the zero width limit. 



> nd{q -my^) , (B.5) 



(g2 - 77^2^)2 + rn^r 



2 r2 



w 



we do obtain that the mean multiplicity is equal to nh{m'^). 

The numerical calculations with the use of formulae (IB. 21) , (IB. 31) result in 
the following values: 

(nw) = 19.04 , (B.6) 

{nt,) = 16.37 . (B.7) 

Thus, the account of non-zero width of the W boson slightly changes the 
average value of the multiplicities. Namely, nw dSlD has gone down by 0.3, 
while ritb ( 155|) has gone up by 0.23, but their sum remains almost unchanged. 
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